Matter-wave solitons in radially periodic potentials 
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We investigate two-dimensional (2D) states of Bose-Einstein condensates (BEC) with self- 
attraction or self-repulsion, trapped in an axially symmetric optical-lattice potential periodic along 
the radius. Unlike previously studied 2D models with Bessel lattices, no localized states exist in the 
linear limit of the present model, hence all localized states are truly nonlinear ones. We consider 
the states trapped in the central potential well, and in remote circular troughs. In both cases, a 
new species, in the form of radial gap solitons, are found in the repulsive model (the gap soliton 
trapped in a circular trough may additionally support stable dark-soliton pairs). In remote troughs, 
stable localized states may assume a ring-like shape, or shrink into strongly localized solitons. The 
existence of stable annular states, both azimuthally uniform and weakly modulated ones, is corrobo- 
rated by simulations of the corresponding Gross-Pitaevskii equation. Dynamics of strongly localized 
solitons circulating in the troughs is also studied. While the solitons with sufflciently small veloci- 
ties are stable, fast solitons gradually decay, due to the leakage of matter into the adjacent trough 
under the action of the centrifugal force. Collisions between solitons are investigated too. Head-on 
collisions of in-phase solitons lead to the collapse; 7r-out of phase solitons bounce many times, but 
eventually merge into a single soliton without collapsing. 

PACS numbers: 03.75.Kk, 42.65.Tg, 42.65.Jx 



I. INTRODUCTION 

Matter-wave solitons, that can be created in Bose- 
Einstein condensates (BECs), are a research subject of 
great interest, both as nonlinear collective excitations 
in macroscopic quantum matter, and due to the poten- 
tial they offer in applications to high-precision interfer- 
ometry, atomic-wave soliton lasers, quantum informa- 
tion processing, and other emerging technologies. Dark 
and bright matter- wave solitons were experimentally cre- 
ated in nearly one-dimensional (ID) trap configurations 
filled with ^'^Rb and ''Li condensate, respectively. 
In these settings, strong radial confinement freezes trans- 
verse dynamics of the condensate, keeping atoms (with 
mass m) in the ground state of the corresponding 2D 
harmonic-oscillator potential, muj'^r'^ /2, while a weak ax- 
ial parabolic trap, allows quasi-free motion of 
the soliton in the axial direction, x. New recent exper- 
iments with ^^Rb 0| and related numerical simulations 
have showed the existence of stable solitary waves in a 
weakly elongated {u!±/uJx — 2.5) trap. These essentially 
three-dimensional (3D) bright solitons exhibit rich and 
complex behavior, not present in their nearly- ID coun- 
terparts. Notable among the new features, are inelastic 
collisions between solitons, which are strongly sensitive to 
phase relations between them and their relative speed. 

Adding a periodic optical-lattice (OL) potential in the 
axial direction of the quasi- ID trap makes it possible to 
create bright matter-wave solitons of the gap type in re- 
pulsive condensates 0, which is a particular manifes- 
tation of very rich physics of BEC trapped in OLs, as 



reviewed in Refs. 0- 

Recently, 2D and 3D localized states in the Gross- 
Pitaevskii equation (GPE, which is the basic mean-field 
BEC model), that may be sustained and stabilized by 
cellular potentials, i.e., multidimensional OLs, have at- 
tracted a great deal of interest, see review A great 
challenge to the experiment is creation of multidimen- 
sional matter- wave solitons, as well as the making of spa- 
tiotemporal solitons in nonlinear optics 8]. In particular, 
multi-dimensional solitons trapped in a low-dimensional 
OL (i.e., ID lattice in the 2D space 0, and 2D lattice in 
the 3D space 0,0|), have been predicted 0. As these 
solitons keep their mobility in the free direction, the lat- 
ter settings can be used to test head-on and tangential 
collisions between solitons. Note that the quasi-lD lat- 
tice potential cannot stabilize 3D solitons 9J (this be- 
comes possible if the quasi- ID lattice is combined with 
the periodic time modulation of the nonlinearity provided 
by the Feshbach-resonance-management technique [Tl| ; a 
general account of the technique of periodic management 
for solitons can be found in book |1^). 

The low-dimensional OL can also be implemented as a 
radial (axially symmetric) lattice in both 2D and 
3D (21| settings. In both cases, stable solitons in the self- 
attractive medium were predicted, in the form of a spot 
trapped either at the central potential well, or (in the 
2D case) in a radial potential trough (additionally, the 
so-called azimuthons were predicted as azimuthally peri- 
odic deformations of vortices in the uniform 2D medium 
^2]; however, they are unstable in the case of the cubic 
nonlinearity). In the latter case, the spot-shaped soli- 
ton may run at a constant angular velocity in the trough 
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(thus being a rotary soliton The circular motion 

may be more convenient for the experimental study of 
mobility and collisions of matter-wave solitons (and soli- 
ton trains) than the already realized soliton-supporting 
settings in the cigar-shaped traps H,!^,^], as in the circu- 
lar geometry the motion is not affected by the weak lon- 
gitudinal trap. This setting is also convenient for studies 
of patterns in self-repulsive BECs. In particular, dark 
solitons were, thus far, created in cigar-shaped traps 0, 
where the background density is essentially non-uniform, 
as it vanishes at edges of the trap. In a circular trough, 
pairs of dark solitons can be created without this compli- 
cation. In fact, both axisymmetric (azimuthally uniform) 
vortices Q and dipole and quadrupole states, that may 
be regarded as stable complexes of two or four dark soli- 
tons |l5i |. were predicted in the two-dimensional GPE 
with the repulsive nonlinearity and Bessel-lattice radial 
potential. The stable dipole and quadrupole patterns 
may also rotate at a constant angular velocity . 

Radial potentials in the form of — eJn (kr) or 
-sJf{kr), adopted in Refs. and [11 UM, respec- 
tively (e and k are positive constants, and Jq and Ji 
are the Bessel functions), can support patterns of most 
types considered in those works in the linear limit of the 
model. Indeed, axisymmetric states trapped in the cen- 
tral potential well, as well as azimuthally uniform vor- 
tices and dipolar and quadrupolar modes trapped in cir- 
cular potential troughs, may be found as 2D eigenstates 
of the quantum-mechanical linear Schrodinger equation 
in such radial potentials. Accordingly, nonlinear states 
with the same symmetries found in Refs. 0, [T3. ITsI l 
may be considered as continuations of their linear coun- 
terparts. An exception is the spot-shaped soliton in a 
circular trough (the one that can perform rotary motion 
along the trough, as mentioned above). On the contrary 
to that, in this paper we report results for the radially pe- 
riodic OL potential, — ecos(2/er), which, obviously, can- 
not support any radially localized state in the linear limit. 
Therefore, all localized patterns reported in this paper 
are essentially nonlinear objects, and they may be differ- 
ent from their counterparts found in the Bessel lattices, 
as concerns their existence and stability conditions. The 
solitons that we report in this work may be categorized 
as solitons of the ordinary type in the model with self- 
attraction, and radial gap solitons (a new species of soli- 
tary waves), in the case of self-repulsion. Note that soli- 
tons of the gap type cannot exist in the above-mentioned 
Bessel-lattice potentials. 

As for the attractive model, a somewhat similar prob- 
lem was considered in Refs. I(i and 0|, where it was 
demonstrated that the self-attractive nonlinearity in the 
ID and 2D geometry (the latter setting is axisymmetric) 
can turn a quasi-bound state in a potential well dug on 
top of a potential hill into a true bound state (including 
2D vortices ^3)- the 3D setting, the situation is more 
complicated, but a similar effect is also possible [l6). 

It is relevant to mention that, besides the attractive 
or repulsive cubic nonlinearity, the radial-lattice poten- 



tial can also be combined with the self-focusing sat- 
urable nonlinearity, characteristic of photorefractive me- 
dia. It was predicte d uM and very recently experimen- 
tally demonstrated flQl that stable solitons are possible 
in this case too (in fact, Ref. 19] presents the first 
ever experimental realization of solitons in any radial 
lattice). Another very recent experimental result is the 
observation of light localization, in the form of spot- 
like and ring-like patterns, in a photorefractive medium 
equipped with a third-order Bessel lattice subjected to 
aziniuthal modulation, which corresponds to a combina- 
tion of saturable nonlinearity and an effective potential 
~ J3 (r/ro)cos'^ (30), where 9 is the angular coordinate 

A different experimental setup that also leads to a 
ring-shaped trap for the BEG has been recently reported, 
based on a specially designed configuration of the exter- 
nal magnetic field that creates a guide for matter waves 
in the form of a torus j2||. In parallel, nearly-lD solitons 
in toroidal traps loaded with self-attractive BEG were 
studied theoretically [23] (in the experiment of Ref. , 
the condensate was self-repulsive). 

The objective of this paper is to investigate various 
types of matter-wave solitons supported by the periodic 
radial lattice potential. It should be mentioned that, be- 
sides the BEG context, the model and results obtained in 
it may also be relevant to photonic-crystal fibers with a 
concentric (rather than usual hexagonal) structure, sim- 
ilar to that in fibers with the multilayer cladding [2J| . In 
that case, the soliton will be a self-trapped beam propa- 
gating in the fiber. 

The paper is organized as follows. In Section II 
we formulate the model based on the respective two- 
dimensional GPE. Section III is dealing with solitons 
trapped in the center (including the radial gap solitons, 
in the case of the model with repulsion), for which re- 
sults are obtained by means of the variational approxi- 
mation (VA; both static and dynamical variants of the 
VA are presented), and in direct simulations. Section 
IV treats solitons trapped in annular channels (potential 
troughs). By means of the VA, we derive an effective 
ID equation for the attractive model, which predicts a 
possibility of the existence of stable axisymmetric ring- 
shaped states, and the onset of their instability against 
azimuthal modulations at a finite density. Exact mod- 
ulated (cnoidal-wave) solutions, which appear above the 
instability threshold of the uniform state, are found too, 
as well as exact solutions for azimuthal solitons. Direct 
numerical simulations corroborate the existence of stable 
uniform and modulated ring-shaped patterns, and reveal 
the radial gap solitons whose central love is trapped in the 
annular trough (in the repulsive version of the model). 
Also reported are stable ring-shaped gap solitons carry- 
ing a pair of dark solitons on their crests. Section V 
addresses dynamics of moving two-dimensional solitons, 
trapped in the circular troughs of the attractive model. 
It is shown that the slow solitons are stable, while fast 
ones are destroyed due to leakage of matter into the outer 
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trough, under the action of the centrifugal force. Head- 
on coUisions between in-phase sohtons lead to collapse 
(intrinsic blow-up), while 7r-out-of-phase solitons collide 
many times, and eventually merge into a single soliton 
(without collapsing); the same happens with in-phase 
solitons colliding tangentially, as they are trapped in ad- 
jacent troughs. Finally, Section VI concludes the paper. 



e ^^*U{r), with a real function U obeying the equation 

^^u+(^ + -^^u-v{T)u + xU^ = Q■ (3) 

\ d?'^ T dr J 

For given potential Vir)^ solutions to Eq. © form fam- 
ilies parameterized by \x or, alternatively, by the norm. 



II. THE MODEL 



The starting point is the standard GPE in a normalized 
form mil, 



du 



ld_ 



1 ^2 \ 

1 u — V{r)u + xl^pu = 0, 



(1) 



written for the single-atom wave function u in polar co- 
ordinates ir,0). Here, X — +1 ^^'^ ~1 correspond to 
the self-focusing and defocusing nonlinearity, respectively 
(alias negative and positive scattering length of inter- 
atomic collisions), and, as said above, V{r) = ecos(2fcr) 
is the radially periodic potential with strength s and 
wavenumber 2k (by dint of obvious rescaling, we set 
k = 1). In terms of OLs, this potential can be created 
by a cylindrical beam whose amplitude is modulated as 
cos(A;r) (see further discussion below). Note that the po- 
tential with e > gives rise to a local minimum at r = 0, 
see an example in Fig. ^ 

Together with Eq. we will use its variational rep- 
resentation. Obviously, the equation can be derived from 
the Lagrangian 



L = 



rdr 



2Tr 



dd 



' du , du* 



du 


2 


^ 1 


du 


dr 


H 




de 



V{r)\i 



,(2) 



where * stands for the complex conjugate. Stationary 




FIG. 1: An example of the radially periodic axisymmetric 
potential, V{r) = — cos(2r). 

axially symmetric states with chemical potential fi are 
looked for as solutions to Eq. ^ in the form of u = 



N ^2tt iU{r)y rdr = N^2TT / \u\^rdr. (4) 
Jo Jo 



Solutions to stationary equation ^ with X ^ +1 
(self-attraction) were obtained by means of a recently 
developed spectral-renormalization method for finding 
self-localized states of nonlinear-Schrodinger (NLS) type 
equations |26j . Independently, the same stationary so- 
lutions were also generated by dint of the well-known 
method which uses the integration of GPE in imag- 
inary time . Stability of the solitary waves was then 
tested by adding small perturbations to them and inte- 
grating the GPE in real time. A stable perturbed soliton 
would shed off some radiation, which was absorbed at 
boundaries of the integration domain, and relax into a 
slightly different stationary form, corresponding to the 
smaller norm. 

In the case of repulsion (x = —1), self-localized states 
were obtained by solving Eq. in real time, with an ab- 
sorber placed at the boundary of the integration domain. 
In that localized waveform (e.g., Gaussian) with 

a suitable norm was chosen as the initial condition. In 
the course of the evolution, excess norm is radiated away 
with linear waves, which are absorbed at the boundary, 
and a self-localized state emerges (see an example in Fig. 
0] below). It is obvious that solitons found this way are 
stable. 

The above-mentioned radial lattices of the Bessel type 
can be naturally generated by a nondiffracting linear op- 
tical beam with the cylindrical symmetry (Bessel beams 
themselves are usually generated by means of axicoms) 
j29l |. On the other hand, the radial potential of the form 
cos(2fcr) can be induced, as said above, by a cylindri- 
cal beam whose amplitude is modulated as cos(fcr). The 
modulation can be provided by passing the cylindrical 
laser beam through a properly shaped plate T^, . Un- 
like the Bessel beam, the one with the cos(fcr) trans- 
verse modulation will suffer conical diffraction, but it is 
not a problem for BEG experiments, as a tight optical 
trap created in the transverse direction may readily con- 
tain a 2D pancake-shaped configuration of the conden- 
sate, of the thickness ~ 2 /im 30]. The diffraction of 
the paraxial beam, with waist diameter D ^ 100 /im, 
relevant to the experiment, on such a short propaga- 
tion distance is completely negligible, as the respective 
diffraction length (alias Rayleigh range) is estimated to 
be Zdiffr / X ^ 10 mm, where A ~ 1 /im is the beam's 
carrier wavelength. 
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III. SOLITONS TRAPPED AT THE CENTER 



Numerical solutions 



It is known that a radial potential structure can easily 
trap a soliton in the central potential well . We start 
the consideration of the present model with solutions of 
this type. Typical examples of the solitonic shape in the 
model with self-attraction (x = +1) are displayed in Fig. 
|5J These solutions were obtained by means of both the 
above-mentioned spectral-renormalization method 
and integration in imaginary time |27| . In the repulsive 
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FIG. 2: Axial cross sections of two stable solitons with equal 
chemical potentials, fi = —0.5, trapped (in the model with 
self-attraction) at the center of the radial structure, V{r) = 
£cos(2r), with e = 1.5 (in this case, the norm is A'^ = 8.9; 
solid line), and e = 0.5 (A'' = 10.1, dashed line). 



model (x — —1), numerical solutions were obtained, as 
said above, by solving Eq. |^ in real time, in the pres- 
ence of the boundary absorber. An example of the relax- 
ation of an initial Gaussian pulse into a soliton is given in 
Fig. 121 Full views of the typical stable solitons trapped 
at the center of the radial structure in the model with 
self- attraction and self-repulsion are displayed in Fig. ^ 

As concerns the solitons in the repulsive model (x — 
— 1), a typical example of which is displayed in the right 
panel of Fig. 21 they may be clearly classified as solitons 
of the gap type, since they exist solely due to the interplay 
of the lattice potential and self-defocusing nonlinearity. 
The ordinary setting that gives rise to gap solitons is 
quite similar, with the OL in the ID, 2D or 3D Cartesian 
coordinates Slowly decaying fringes, attached to the 
central core of the soliton and evident in the right panel 
of Fig. 01(and absent in the ordinary soliton displayed in 
the left panel), is a characteristic feature of gap solitons. 
Thus, this self-trapped state may be naturally called a 
radial gap soliton. 
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FIG. 3: Formation of a self-trapped state {radial gap soliton) 
in repulsive BEC in the radial lattice with strength e — 3 (in 
the presence of absorbers at the boundary of the integration 
domain) from an initial Gaussian with norm N = Gn. The 
norm of the eventually established state is A' = 5.63, i.e., 
« 30% of the initial value. 




FIG. 4: Left panel: The same localized state in the attractive 
BEC whose axial cross section is shown by the solid curve in 
Fig. H(with N = 8.9, ^ = -0.5, and e = 1.5). Right panel: A 
localized state in the repulsive BEC, with A'^ = 5.63, fi = 1.88, 
amplitude A — 1.358, and e = 3, which was generated as 
shown in Fig. |21 



B. Variational approximation 

To describe a stationary soliton trapped in the central 
potential well in an analytical approximation, we adopt 
a Gaussian ansatz. 



U{r) — A exp 



r 

2^ 



(5) 



Comparison with Fig. (21 shows that, in the model with 
attraction (x = +1), this ansatz is quite adequate for 
solitons with a larger norm, trapped in a relatively weak 
lattice, but it may be inaccurate as an approximation for 
solitons with smaller TV trapped in a stronger lattice. 

The application of the variational approximation (VA) 
to stationary equation (j^J, with the help of the accord- 
ingly modified Lagrangian is straightforward (cf. 
Refs. 0,112, 1131, where the same Gaussian ansatz was 
used to predict 2D solitons trapped in quasi-lD and 
square OLs). Thus we arrive at a set of VA-generated 
equations that relate the soliton's size a to its norm and 
chemical potential: 



47r 



1 



1 - odf 



da 



-4' (6) 
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where /(a) = \ — ^ae~°' erfi(a), erfi(a) = erf(za)/i, and 
erf is the standard error function. 

The analysis of dependence = ^^{N) following from 
Eq. © with X = +1 (self- focusing) predicts the exis- 
tence of a family of solitons trapped at the center of the 
radial lattice, which are stable accordin g to the Vakhitov- 
Kolokolov (VK) stability criterion [3lT36||. dn/dN < 0, 
as shown in Fig. |S1 Although the VK criterion is only 
necessary for the stability, as it ignores a possibility of os- 
cillatory instability with complex eigenvalues, direct sim- 
ulations of Eq. confirm the stability of these solutions 
(see below). 




FIG. 5: Dependences ^ = /i(A^) following from variational 
equations (|SJ, for the self-attraction case, x = +1- According 
to the VK criterion, dfi/dN < 0, stable solutions belong to 
upper branches of the curves. 

Additional evidence for the existence of self-trapped 
localized states of BEC in the present model can be pro- 
vided by a more general, time-dependent, version of the 
VA. To this end, we define a generalized ansatz [cf. Eq. 
©] 

u{r, t) — A[t) exp 



2{a{t)y 



(7) 



where 6 is a real radial chirp. Applying the standard 
VA procedure Isij , one arrives at the following evolution 
equation for the width of the localized state, 

d^a 4(1 -x') 



dt^ 



+ 2e 



2a + V7r(l-2a2)e"'' crfi(c 



(8) 
(9) 



where an effective nonlinearity strength is 

and the amplitude is given by 

A^{t)^N/{Tra^{t)). (10) 

Equation |(SJ describes motion of a unit-mass particle 
with the coordinate a{t) in an effective potential 

2(1 -X') 



In Fig. I^lthe effective potential is depicted for differ- 
ent values of x': for the attractive and repulsive versions 
of the model. 




ID 1 5 2 25 3.0 3 5 40 




1 5 2.Q 2.5 3.0 3.5 



FIG. 6: The variational potential, as given by Eq. Ijll^ for 
different values of x = X-^/(47i"), at fixed strength e of the 
radial optical lattice. 

Solitons in the attractive BEC (x = +1) exist when po- 
tential 1)11(1 possesses a local minimum. In particular, for 
e = 0.5 (see Fig. EJ, the potential minimum exists pro- 
vided that x' > 0.45. For the norm smaller than that cor- 
responding to x' = 0.45 [see Eq. ©J, the self-attraction 
is too weak to form a soliton, and the matter-wave pulse 
spreads out through barriers of the radial potential. On 
the other hand, in the case of attraction, x' is also limited 
from above, by the collapse occurring at a critical value, 
Xcr = ^cr/(47r). In the absence of the lattice (e = 0), 

a numerically exact critical norm is iVi°^ ~ 11.7 (it IS 
equal to the norm of the Townes soliton) |3(i| . For any e, 
the VA predicts Xcr = 1, as the value of coefficient x' at 
which the first term in effective potential ((ll|l vanishes. 



(var) 



= 47r 



According to Eq. ©, this gives N( 

Thus, the VA predicts a finite existence region for sta- 
ble self-trapped localized states in the attractive BEC; 
for example, it is 0.45 < x' < 1 for e — 0.5. While the 



47r, nu- 



Uia) 



2i/7reae " erfi(a). 



(11) 



upper edge of this region is fixed at = Ai™ 
merical analysis of expression (|ll|l demonstrates that the 
lower critical norm necessary for the existence of the self- 
trapped state almost linearly decreases with the increase 
of the OL strength e. Numerically exact solution of the 9- 
independent version of full CPE Q demonstrates that, 
although the collapse sets in at a critical value of the 
norm which is slightly higher than Ai™'-* = 47r, similar 
to the case of the 2D soliton stabilized by the quasi- ID 
or square OL [s^, HI] , the dependence of Ncr on e is 
weak. The physical mechanism behind the rise of the col- 
lapse threshold in a trapping potential is that the confine- 
ment allows the soliton to increase the internal "quantum 
pressure" , generated by the dispersion term in CPE Q , 
which counteracts the nonlinear self-focusing of the wave 
packet. 

Predictions of the VA for dynamical regimes can be 
compared to direct simulations in terms of the time de- 
pendence of the amplitude, A(t), in a soliton which was 
excited by a sudden perturbation [the variational predic- 
tion for A{t) is taken as per Eq. (|10|l ]. To display an 
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example, we take the stationary solution from Fig. 13 
corresponding to e = 0.5 and N — 10.1, and suddenly 
increase the lattice strength to e — 0.75. Comparison 
of the variational results with simulations of Eq. ^ is 
presented in Fig. |7| As observed, agreement for absolute 
values of the amplitude is, at best, qualitative, while the 
oscillation frequencies coincide much better. 



GPE 




VA 



FIG. 7: Comparison of a result of the variational approxima- 
tion ("VA") with direct simulations of Eq. ((TJ ("GPE") with 
X = +1, in terms of the soliton's amplitude as a function of 
time. Intrinsic vibrations of the soliton were triggered by sud- 
den increase of the lattice strength from e = 0.5to£ = 0.75, 
for the self-trapped state with = 10.1 in Fig. |21 

In the case of the repulsive BEC (x = —1), the behav- 
ior is opposite, as seen in the right panel of Fig. In this 
case, a larger norm for given OL strength e causes spread- 
ing of the localized state, because the self-defocusing non- 
linearity overcomes the confining effect of the trapping 
potential. Simulations of Eq. |^ with x — — ^ show that 
the wave packet sheds off excessive norm in the form of 
linear waves, which are absorbed at domain boundaries, 
and relaxes into a stable shape after the norm becomes 
adequate. For example, at e = 2 stable localized states 
in the repulsive BEC exist at \x'\ ^ 1 [i.e., the norm is 
limited to values N < An, see Eq. ©]. 



IV. RING-SHAPED SOLITONS 
A. Analytical approach 

The most essential property of radially periodic poten- 
tials is that localized states can self-trap not only in the 
central potential well, but also in remote circular troughs 
[l3|. If the trough's curvature is not significant, these lo- 
calized states resemble 2D solitons in the quasi-lD OL re- 
ported in Ref. Here, we aim to derive an asymptotic 
ID equation for patterns trapped in a circular trough, 
and obtain its relevant solutions. To this end, we adopt 
the following ansatz for the solution. 



ip{r, e, t) = V2A{y, t)sech (xjp) , 



(12) 



where y = vqO is the coordinate running along the cir- 
cumference of radius rg, A(y,t) is a slowly varying com- 
plex amplitude, and the transverse variable is a; = r— 7rm, 



with To = Trm being a radial potential minimum in a 
vicinity of which the ring-shaped pattern is trapped. As- 
suming m to be large enough, we consider the annular 
pattern placed in a remote circular trough, which may be 
treated as a quasi-rectilinear potential trap, with negli- 
gible curvature. Obviously, the latter condition amounts 
to p <C ro [recall p is the radial size in ansatz (|12|l ] . 

Neglecting the curvature in Eq. JQ, we rewrite it as 
a 2D equation with x and y treated as local Cartesian 
coordinates. 



9y2 



''^ + in^ + l^+x\4'\'i'-ecosi2x)^ = 0. (13) 



The next step is to reduce Eq. H13() to an effectively ID 
equation along azimuthal direction y. To this end, we 
note that Eq. (|13|l corresponds to the Lagrangian [cf. 
Eq. ©] 



L 



dxdy 



(14) 



'^X \ipf - ecos(2a;) 



In this Lagrangian, we substitute ansatz (|12|l and inte- 
grate the resulting expression in transverse direction x. 
After a straightforward calculation, this yields an effec- 
tive (averaged) Lagrangian, 



/ + 00 
dy 
-OO 



ipiAtA* - AA*t) - 2p - 2I^\A\^ 



2Mr 

3p 



:XP\A\ 



4 2Trep^ \A\^ 
sinh (tt/j) 



(15) 



sech'*x) x^dx 



0.61. 



where I = °° (sech' 

The application of the standard variational procedure 
to averaged Lagrangian (|15(l gives rise to a cumbersome 
system of coupled equations for the complex amplitude 
A{y,t) and real width p{y,t). The equations strongly 
simply if we assume that A and p are slowly varying 
functions of y, which, in particular, corresponds to long- 
wave perturbations that account for the onset of the mod- 
ulational instability (MI), see below. In particular, in 
the attractive model (% — +1) one may first adopt the 
lowest-order approximation, completely neglecting the y- 
dependence and assumes time- independent |^| and p, 
which gives rise to a result previously known from the 
application of the VA to the NLS equation in ID with 
potential cos(2fca;) |3^ : 



7rp 



cosh(7rp) 1 
sinh^(7rp) sinh (tt/ci) 



(16) 



This relation is cumbersome too, but straightforward 
analysis demonstrates that it takes a very simple form, 
1^1 = 1/p, when £ is both small and large [in the latter 
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case, scaling \A\ ~ ~ e~^^^ is assumed, to lend the 
last three term in the Lagrangian density in Eq. 115|1 the 
same order of magnitude] . 

To obtain a closed-form Lagrangian for the amplitude 
field A{y, t) in the attractive model, one should express 
p in Eq. IjlSfl in terms of \A\^ using relation This 
will lead, in the general case, to a messy result; however, 
the above-mentioned simplest approximation, p = 1/ 
which replaces Eq. (|15|l . yields a tractable expression. 



^off 




[the last term in the Lagrangian density should be 
dropped if e is small; if e is large, the term is obtained 
within the framework of the above scaling, that assumes 
small p, hence sinh(7rp) « 7rp]. Finally, the action func- 
tional, S — J Lcffdt, as defined by Lagrangian l|17|) . gives 
rise to the Euler-Lagrange equation, SS/5A* = 0, in the 
following form: 



iAt + (2 + I) Ay 



3/ 

T 



A 



— A* 

A* 



2A* 



A* 



A 



2 _ 
A 



(18) 



where A{y,t) = A(y,t)e~*^* [this transformation elimi- 
nates a linear term in the equation generated by the last 
term in Lagrangian (|17() ]. 

Formally setting / = in Eq. H18() recovers an equation 
known as a model for the propagation of surface waves 
on a plasma layer with a sharp boundary [39| . A class of 
nonpolynomial generalized NLS equations similar to Eq. 
(|18|l . and some of their fundamental solutions, such as 
single-soliton ones, were introduced in Ref. (see also 
Ref. [41]). 

Stationary solutions to Eq. p8|l with chemical 
potential p are looked for in the ordinary form, 
A{y,t) — e^*''*i?(y), with B{y) satisfying equation (1 + 
/) [2B" - {B'f/B]+2B^ + pB = 0, where B' = dB/dy. 
This equation has its Hamiltonian, 



For H < and /z < 0, a family of periodic cnoidal-wave 
solutions is obtained in the form 



Scn(y) 



BqBi 



(21) 

where cn(z, q) is the Jacobi's elliptic cosine, with modulus 



'Bq-Bi 

Bn 



IB2 



Bi + IB2I 



< 1, 



(22) 



and B2 < < Bi < Bq are three roots of equation 
(2/3)5^ + pB - H = [it is obtained by setting B' = 
in Eq. ^^]. Solutions (|2U exist, for given p < 0, 
in a region of < —H < (\/2/3) l^l"^/^; as said above, 
the limit oi H ^ corresponds to soliton (|2U|I . and the 
opposite limit, H ~ ~ (•\/2/3) |^|'^/^, corresponds to a 
uniform CW (continuous-wave) solution, with 



Bo = Bi = B{y) = 



(23) 



Cnoidal solutions on the ring of radius vq must satisfy 
the corresponding periodic boundary conditions (b.c), 
A{y,t) = A{y + 27rro,t) (for the ordinary ID NLS equa- 
tions with repulsion and attraction, analytical solutions 
satisfying the periodic b.c. are collected in Refs. 
and |43j| . respectively). Accordingly, parameter H in so- 
lutions (|21|) . 1)22(1 is not continuous, but rather takes dis- 
crete values selected by matching the b.c. to the period- 
icity of cn^. 



vrro 



' 6(1 + /) 
Bo (S1 + IB2I) 



(24) 



where K is the complete elliptic integral, and n is an 
arbitrary integer. 

Note that Eq. H18|l features the Galilean invariance: if 
A{y, t) is a solution, then its counterpart in the form of 
a solution moving at an arbitrary velocity v is 



Ac{y,t) = A{y - ut)exp 



(25) 



hence one can boost solitons ()20|l and cnoidal waves H21|l 
to a (formally) arbitrary velocity by means of transfor- 
mation H25() . In fact, the velocity is not arbitrary because 
the factor exp [ivy /A) also must satisfy the periodic b.c, 
which leads to a quantization condition, v = 4n/ro, with 
n = 0,±1,±2, ... . 



H = (2/3)^3 + pB + {l + I){B'f/B, 



(19) 



Setting H — 0, one obtains a family of soliton solutions, 
for any p < 0: 



Bsoiiy) = 



-— sech 
2 



(20) 



B. Modulational instability of the uniform ring 
soliton in the model with attraction 



In the model with attraction (% — +1), the existence of 
uniform ring-shaped solitons trapped in annular troughs 
of large radii is obvious, while a nontrivial issue is a pos- 
sibility of finding a stability region for such solutions. In 
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terms of Eq. 1)18(1 . the ring soliton is represented by a 
CW (continuous-wave) solution with a constant amph- 
tude Bq, which was mentioned above as a hmit case of 
the cnoidal solution family, with = —2Bq. 

It is well known that, while all CW solutions of NLS- 
like equations with attraction are modulationally unsta- 
ble [including nonpolynomial equations similar to Eq. 
(|18|l 40]], periodic b.c. may stabilize them if the CW 
amplitude is smaller than a certain critical value. This 
argument suggests a possibility to find a stability region 
for the uniform ring solitons in the present model with 
attraction. 

A standard analysis of the MI assumes a perturbation 
of the CW solution, 

i(y, t) = [Bo + Bi(y, t)] exp {2tB^t + ixi {y, t)) , 

with infinitesimal amplitude and phase perturbations Bi 
and xi. Eigenmodes of the perturbations may be looked 
for assuming that they are proportional to exp (ipy + jt), 
where p is an arbitrary real wavenumber, and 7 is the 
corresponding instability growth rate. Straightforward 
calculations demonstrate that the CW solution to Eq. 
()18|l is subject to the MI under the condition 

p2 [2i?2 _ (1 + /)p2] > 0. (26) 

In the annular system that we are dealing with, p is 
subjected to the geometric quantization, the same way as 
above, i.e., p — n/rQ with integer n. Therefore, condition 
()26|l with lowest \n\ = 1 is not satisfied, making the CW 
solution modulationally stable, under the condition 

2B^rl <l + I. (27) 

In a final form, this stability condition may be expressed 
in terms of the full norm of the axisymmetric ring soli- 
ton. Substituting ansatz H12|l in Eq. Q, and mak- 
ing use of the assumptions adopted above, one obtains 
N w SnropB^. Further, substituting here the above ap- 
proximation, p = 1/Bq, yields 

N = SnroBa. (28) 

With regard to this relation, stability condition H27|) 
amounts to a limitation on N: 

N < iVthr = 4v/2(l + J)7r w 22.55. (29) 

(note that this result does not contain any parameter). 

A natural way to look at these stability conditions is 
to consider a situation with gradually increasing N (or 
-Bo)- The azimuthally uniform ring solution will loose its 
stability when the amplitude (or norm) attains the criti- 
cal value defined by Eq. [or by Eq. . One may 
expect that the onset of the MI gives rise to a bifurca- 
tion, which creates stable azimuthally modulated solu- 
tions, with a modulation depth scaling as ^/N — iVthr for 
< {N — Nthr) /Nthi ^ 1. Alternatively, one may fix the 
CW amplitude, Bq, and gradually increase the trough's 



radius rg ; the loss of the CW stability and emergence of 
the modulated solutions should occur when vq attains a 
critical value following from Eq. (|27|) . 

{ro),,. = Vi^ + I)/2Bo'. (30) 

In fact, the modulated states generated by the MI on- 
set should be nothing else but the stationary cnoidal wave 
(EU, Indeed, with regard to Eq. and the fact 

that K{Q) — tt/2, it follows form matching condition H24() 
that, with the increase of ro for fixed Bq, the cnoidal- 
wave solution appears, with an infinitely small modula- 
tion depth, at rg = (^o);.!-, where the critical value iro)^^. 
is precisely the same as defined by the onset of the MI 
[i.e., given by Eq. (|Sn)l ]. 

A defect of the above approximate analysis is that, 
with regard to relation Aq = l/p, stability condition 
may be cast in the form of p/rg > y/2/{l + I) w 1.1, 
which does not comply with the underlying low-curvature 
assumption, po <^ rp. This fact makes the existence of 
modulationally stable ring-shaped solitons in the present 
model with attraction doubtful. In direct simulations of 
Eq. ^ with X = +!> we were unable to find modula- 
tionally stable annular uniform nonlinear states. In this 
connection, it is relevant to recall that the linear limit of 
Eq. with potential V{r) = ecos(2r) does not allow 
any stationary solution localized in the radial direction, 
therefore the weak nonlinearity is essential in stabilizing 
the axisymmetric rings. 

Note that localized ring-shaped states may be possible 
in the linear limit of the Bessel-lattice models consid- 
ered in Refs. j3 s-^d 0], with potentials — e Jq (fcr) 
and — eJ^(fcr), respectively; however, the azimuthal 
(in)stability of such states was not studied in the presence 
of attraction. On the other hand, stable annular patterns 
(both static ones and persistent ring-shaped breathers) 
were recently found in an attractive model with a sin- 
gle circular trough, which was created not by the linear 
potential, but rather by modulation of the nonlinear coef- 
ficient in the radial direction (4J] . In that model, annular 
states may only be stable if the annulus is wide enough, 
with the ratio of its outer and inner radii exceeding 0.47. 
Vorticity-carrying rings were found to be unstable in the 
same setting. 

C. Gap Soliton nature of ring shaped states 

In this section we investigate the gap soliton nature 
and the band structure associated with ring shaped states 
in terms of the following 2D CPE 

iipt = ~Ari^ + scos{2r)ip + x\M'^ip (31) 

where A^. denotes the spherical 2D laplacian: — 
-I- — ^, with l'^ = — the squared 2D angular 
momentum. This equation corresponds to the ^ ^ f 
limit {u{r, n 12, if) -(/'(r, if)) of the spherically symmet- 
ric 3D GPE, with (r, 6, ip) denoting spherical coordinates. 
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Notice that apart for the factor of two in front of the term 
due to spherical rather than cyhndrical symmetry, 
Eq. H^Ufl is the same as Eq. From a physical point 
of view Eq. (|31|) is appropriate to describe a condensate 
trapped in the outer part of the conical region created by 
a cylindrical beam of light passed through a very short 
focal distance lens. In the following we shall refer to 
this setting as the "pizza" geometry since, in contrast 
with the pancake geometry, the thickness of the conden- 
sate in the axial direction is not uniform and reduces to 
zero as r — > 0. This geometry could be implemented 
by combining a cylindrical beam creating the radial lat- 
tice with a blue-detuned conical beam (the pizza-shaper) 
launched in the opposite direction. Notice that at large 



iij 




n n 



FIG. 8: Left Panel: Energy levels of Eq. ^ with x = -1, 
I = 10 and e — 10. Right panel: The same as in the left panel 
for the repulsive case x = 1 with I = 2 and e — —4. 

radial distances (i.e. in the outer part of the pizza) the 
differences between pancake and pizza geometry become 
negligible especially if the thickness of the border is very 
small compared to the radius. Ring shaped states of Eq. 

located in potential wells far from the center are ex- 
pected, therefore, to be qualitatively similar to the ones 
of the pancake setting Eq. QJ. 

In the following we look for ring gap solitons of Eq. 

in the form of stationary states ip{r, ip, t) of the form 

^(r,^,t) = (/)(r)e-'^*e*''^. (32) 

By substituting Eq. (|32|l into Ea. (|31|l we obtain the fol- 
lowing equation for the radial function 4>{r) 

2 P 
(j)„ + + (^ - £ cos(2r))0 - (x|0p + = 0, (33) 

where I are integer numbers which allow to satisfy the 
periodic boundary condition in ip: 4^{r,ip + 27r,t) = 
-(/^(r, if, t). Equation H33|l can be mapped into a ID eigen- 
value problem by means of the transformation 4>{r) = 
leading to the following reduced radial equation 

Crr + [E-e cos(2r) ^^^§^]^ = (34) 

which has the form of a ID Schrodinger equation with 
periodic potential and nonlinear effective centrifugal bar- 
rier: {P + xl^P)/^^ (here E = For sufficiently large 



distances r ^ 1 and fixed angular momentum and non- 
linearity, this equation reduces to the Mathieu equation 
with the well known band gap structure. The part of 
the spectrum associated to excitations localized nearby 
the origin, however, is expected to be strongly influenced 
by the centrifugal barrier which at short distances en- 
hances the nonlinearity and favors the formation of lo- 
calized states. 

To investigate the band structure and the gap solitons 
of Eq. lIHUl we use a self-consistent method developed 
in Ref. |33- The gap solitons obtained with the self- 
consistent method will then be used as initial conditions 
for time propagation of the radial 2D GPE in Eq. 



5- 




5 10 15 20 25 30 a 5 10 15 2 25 30 



FIG. 9: Panels a,b. Attractive radial gap solitons (thick lines) 
trapped in circular troughs of the potential at r = I.Stt and 
corresponding to the energy levels a, b of the left panel of 
Fig. |H] respectively. Panels c,d. Repulsive radial gap solitons 
(thick lines) trapped in circular troughs of the potential at 
r = I.Stt and corresponding to the energy levels c,d of the 
right panel of Fig. |H| respectively. The thin lines corresponds 
to the radial potential ecos(2r) scaled by a factor 50 in the 
left panel and by a factor 20 in the right panel, for graphical 
convenience. 

As a result we show the existence of radially symmetric 
gap solitons which correspond to the ring shaped states 
discussed before. In Fig. |Slwe depict the lower energy 
part of the spectrum obtained from Eq. I|34|) in the case 
of attractive and repulsive interactions. We see in both 
cases the existence of a series of discrete levels located in 
the gaps between the bands which correspond to radially 
symmetric localized states. The wave functions of the 
energy levels labeled a-d in Fig. (Sjare shown in Fig. |2| 

The symmetries of these localized states is similar to 
the ones of intrinsic localized modes of nonlinear lattices 
and gap solitons of continuous ID GPE with a periodic 
potential. Adopting the same terminology introduced for 
these cases we shall refer to them as the onsite symmet- 
ric and the onsite asymmetric ring gap solitons. A 3D 
plot of these states is reported in Fig. ^| (notice that 
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FIG. 10: Top panels. 3D view of the gap solitons correspond- 
ing to the bound states a (left panel) and b (right panel) of 
Fig. |H] Bottom panels. 3D view of the gap solitons corre- 
sponding to the bound states c (left panel) and d (right panel) 

of Fig. li 




FIG. 11: Time evolution {y = 0) section of \-ip\'^ obtained from 
Eq. Q by taking as initial conditions the localized states de- 
picted in panels a-d of Fig. El The norm of the attractive 
states in panels a,b are 4.85, 4.81, respectively, and corre- 
spond to bound states a,b, of Fig. |H| Repulsive states in 
panels c, d have norm 11.47 and 9.17, respectively, and corre- 
spond to the bound states c,d, of Fig. |H| Parameters are the 
same as in corresponding panels of Fig. |S] 

the asymmetric state gives two rings of matter located 
in the same trough of the potential in the attractive case 
and many ring oscillations in the repulsive case ). The 
stability of these gap solitons has been checked by time 
propagation under the 2D GPE in Eq. The re- 

sults are reported in Fig. for the y — section of 
of \ip{x,y,t)\'^ as a function of time. We see that for 
the attractive case the gap solitons are always unstable, 
the onsite symmetric state becomes modulationally un- 



stable at t « 80 while the asymmetric state decays into 
a state localized in the center and a series of ring soli- 
tons which expand trough the lattice keeping the same 
symmetry as the initial state. For the repulsive case we 
see that while the asymmetric state in the second band 
gap (panel d) is unstable, the onsite symmetric ring gap 
soliton in the first band gap is quite stable under time 
evolution. The stability of the onsite symmetric repulsive 
ring gap solitons and the modulational instability of the 
repulsive ones has been found also for other parameter 
values and seems to be a general property of the model. 
The existence of other types of ring gap solitons (which 
resemble intersite symmetric and asymmetric intrinsic lo- 
calized modes of nonlinear lattices) is also possible but 
they look even more unstable under time evolution than 
the onsite asymmetric ones for both attractive and repul- 
sive cases. In the following section we shall investigate 
onsite symmetric gap soliton states in more details. 



D. Numerical investigation of ring-shaped states 

Generic examples of a ring-shaped soliton in the at- 
tractive model, and its counterpart in the repulsive one, 
found from direct solution of GPE in, respectively, 
imaginary and real time, are displayed in Fig. 1121 (in 
the latter case, real-time simulations were run with ab- 
sorbers set at r = Stt and r = 0). The latter example 
represents radial gap solitons with the core trapped in 
a remote circular trough (cf. the right panel in Fig. ^ 
that shows a radial gap soliton with the core trapped in 
the center). The fact that the soliton belongs to the gap 
type is clearly attested to by conspicuous inner and outer 
radial fringes attached to the core (the fringes are absent 
in the ordinary soliton in the left panel of the figure). 




FIG. 12: Left panel: An unstable ring-shaped soliton in the 
attractive model with £ = 2, trapped in a circular trough at 
ro = 67r. The amplitude of this soliton is 1.8, and its chemical 
potential and norm are fj, — —2.3 and A'' ~ 420. Note that 
analytical approximation l|28^ yields, in this case, N ~ 426, 
which is very close to the numerically found value, and its 
modulational instability is definitely predicted by Eq. (1291 . 
Right panel: A stable ring-shaped soliton (actually, it is a 
radial gap soliton) in the repulsive model with e = 3, trapped 
at ro — 5n. In the latter case, A'^ — 344.3 and /j, — 0.78. 

The evolution of the unstable axisymmetric ring from 
Fig. El (left panel), and outcome of the evolution are 
displayed in Fig. E| It is observed that the system tends 
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to form a stationary necklace-like pattern composed of 
six peaks towering above the remaining quasi-uniform 
background. 




FIG. 13: (a) Development of the modulational instabil- 
ity of the axisymmetric state in Fig. 1121 (left panel), 
triggered by a small azimuthal perturbation, 5u = 
0.02 cos(6(^) exp (— (r — ro)^/2). The evolution is shown in 
terms of the amplitude, taken along circumference 

y — roO, with ro = 6n, 6 running from to vr (due to the sym- 
metry, only half of the circle is shown), (b) The field profile 
formed by t = 5, demonstrating the emergence of a regular 
necklace-shaped pattern. 

In the repulsive model, a stable ring (i.e., the radial 
gap soliton) may additionally carry, on its crests, pairs 
of dark solitons. An example of such a pattern is dis- 
played in Fig. (it resembles the dipole-mode ring soli- 
ton recently found in the model with repulsion and Bessel 
lattice in Ref. 15]). 

Stable static rings with a density modulation in the az- 
imuthal direction, corresponding to analytical solutions 
(|21|l for the cnoidal waves, have been found too in the 
attractive model. A typical example of such a pattern is 
displayed in Fig. in fact, it was generated as an out- 
come of the evolution of an unstable axisymmetric ring 
with the same norm. 

Increase of the norm leads to a transition from weakly 




FIG. 14: An example of a stable configuration formed by two 
dark solitons created on top of a ring-shaped soliton in the 
repulsive model with e = 3, trapped at ro = Stt. This state 
has norm A'^ = 330.6 and chemical potential ^ = 0.80. 




FIG. 15: A stable ring-shaped soliton featuring weak az- 
imuthal modulation, in the model with attraction and strong 
radial lattice (e = 10). The solution's norm and chemical 
potential are = 7.2 and /i = —5.7. 

modulated stable patterns to deeply modulated ones. An 
example of such a stable configuration, generated from 
four Gaussian pulses, with the phase shifts tt between 
them, by integration of Eq. ^ in imaginary time, is 
displayed in Fig. ^| (note that the total norm of this 
pattern exceeds that corresponding to Fig. Elby a factor 
of ~ 3.5). 

V. ROTATIONAL DYNAMICS AND 
COLLISIONS BETWEEN SOLITONS 

A. Stability of rotary solitons 

Strongly localized solitary waves of the attractive BEG, 
self-trapped in a large-radius annular potential channel. 
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FIG. 16: A stable set of four solitons in the attractive model, 
trapped in a circular trough with ro = 5n and e = 2. The 
norm of each soliton is A*' = 27r, with the phase difference tt 
between adjacent ones. The lattice is embedded in an inte- 
gration domain of size IGtt x IGtt. 

like those shown in Fig. 1161 can freely move along in the 
channel. If set in motion with a sufficiently small velocity, 
the soliton can circulate in the channel indefinitely long, 
preserving its integrity, as shown in Fig. 1171 




FIG. 17: A stable soliton in the attractive model, with norm 
— 9.42, trapped in the circular trough at ro — 5tv, was set 
in motion in the clockwise direction with velocity v = 1. As 
illustrated by the sequence of snapshots, the circular motion 
in the trough continues indefinitely long without any visible 
loss. 

However, when the circulation speed is too large, the 
centrifugal force acting on the soliton can cause the mat- 
ter to tunnel into adjacent radial troughs. The loss of 
the norm (number of atoms) caused by the centrifugal 



tunneling eventually brings the soliton's norm below the 
threshold value necessary for its existence, which leads 
to disintegration of the soliton, see Fig. ITSl 




FIG. 18: If the same soliton as in Fig. 1171 (with initial norm 
= 9.42) is set in motion with a greater velocity, n = 3, 
the centrifugal force gives rise to the underbarrier leakage of 
atoms into the outer circular troughs. As illustrated by the 
sequence of snapshots, the decaying soliton cannot preserve its 
shape. The norm remaining within the given radial channel 
at t — 400 is A'^anai — 0.5, i.e., ~ 5% of the initial value. 

If the lattice is weak, and the soliton is strongly self- 
trapped (with a sufficiently large norm), the centrifugal 
force can also drive it away from the center, across the 
lattice. In that case (not shown here), the soliton also 
suffers the radiation loss, and eventually disintegrates. 

B. Collisions between solitons 

The present model offers unique possibilities for ex- 
ploring interactions and collisions between matter-wave 
solitons, as the dynamics in circular channels is free of 
external perturbations, which are always present in pre- 
viously reported settings (in the form of the weak longi- 
tudinal trap intended to keep the solitons within a finite 
spatial domain). The same advantage is offered by se- 
tups with the toroidal trap [2^1 123|. However, sideline 
(tangential) collisions, which are possible between soli- 
tons moving in adjacent potential troughs of the radial 
lattice, have no counterparts in the toroidal traps. 

In head-on collisions of solitons in the attractive model 
(in the same circular channel), two competing processes 
play major roles. The first is the nonlinear self-focusing, 
which may give rise to intrinsic collapse of a large-norm 
"lump" temporarily formed by the colliding solitons, and 
the other is the interference due to coherence of the 
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matter- wave solitons [9|, |33|. The interference pattern, 
featuring alternating regions of high and low density, can 
suppress the collapse. The predominance of either mech- 
anism depends on the time of interaction. If it is small 
(in the case of large collision speed), solitons can pass 
trough each other without collapse (although with con- 
spicuous loss, see below, hence the collision cannot be 
termed quasi-elastic), even if their total norm exceeds 
the collapse threshold. On the contrary, slow collisions 
lead to the onset of collapse. These two possible scenar- 
ios of the head-on collisions between in-phase (mutually 
attracting) solitons are shown in Figs. 1191 andEUl 

Figure El presents an example of a relatively fast col- 
lision. At the moment of the full overlap between the 
solitons, t — 4.3, the interference fringes are evident, and 
are also seen in the cross section along the circumference 
of the channel (lower right panel). Although the total 
norm of the colliding solitons is overcritical (A'tot = 14.9), 
and the collision time {At ~ 0.5) is larger than the col- 
lapse time (which is estimated as icoiiapse ^ l/-^tot) , the 
solitons separate without collapse. Nevertheless, the col- 
lision is essentially inelastic, the loss being additionally 
enhanced by the matter leakage under the action of the 
centrifugal force. On the other hand. Fig. 1201 shows that 
the solitons colliding at a sufficiently small velocity merge 
together and indeed blow up due to the collapse. 
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FIG. 19: Collision between two solitons with individual norms 
= 2.5n and zero phase difference, in the circular trough at 
ro = 5tt. One of the solitons is set in motion (at x = 57r, y = 0) 
at f = with speed n = 3 in the clockwise direction, while 
the second one (at x — 0,y = —5n) stays quiescent. The 
norm remaining within the given circular channel at f = 10 is 
Affinal = 12.75. The soliton which was originally at rest loses 
less norm than the moving one. 



Two TT out-of-phase solitons colliding in the same po- 




FIG. 20: The head-on collision of two in-phase solitons in the 
attractive model, when both of them were set in motion with 
relatively small velocities, v = ±1, in one circular channel 
(ro — 5n), ends up with collapse. Upper panels: density plots 
at t = and the moment of collapse {t = 12). Lower panel: 
the evolution of the density profiles in the cross section along 
the circumference of the potential trough. As in Fig. 1191 the 
initial norm of each soliton is A'' = 2.57r. 



tential trough bounce from each other due to repulsion 
between them. Then, repeated bouncing collisions are 
observed, which is similar to the dynamics in the Bessel 
lattice ^ a quasi-lD OL, equipped with 

a weak longitudinal trap jSj. However, in contrast to 
those works, we have observed that many collisions in the 
circular trough gradually wash out the phase coherence 
between the solitons. As a result, the solitons cease to 
repel each other, and a clear interference pattern disap- 
pears. Eventually, the originally repelling solitons merge 
together, shedding the excessive norm away, as shown in 
Fig. 1211 The intrinsic collapse does not occur in this case, 
as the loss makes the norm of the finally established sin- 
gle soliton smaller than the critical value at the collapse 
threshold. 

Collisions of matter-wave solitons in the attractive 
BEG is accompanied by strong exchange of matter be- 
tween them. The present setting makes it possible to 
investigate related phenomena for solitons moving in ad- 
jacent channels of the radial potential, and thus colliding 
tangentially. The result is that two in-phase solitons are 
unstable against the flow of matter in this setting, and 
tend to merge (a similar trend was observed in the model 
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VI. CONCLUSIONS 



In this work, we have proposed a new setting to explore 
two-dimensional (2D) localized states in self-attractive 
and repulsive BECs, in the form of a periodic radial op- 
tical lattice. A crucial difference from previously studied 
2D models with Bessel lattices is the fact that no local- 
ized state may exist in the linear limit of the present 
model, hence all solitary states in it are "true solitons" , 
impossible without the nonlinearity. Besides the BEC, 
the model may also apply to spatial solitons (beams) in 
photonic crystal fibers with a circular intrinsic structure. 



FIG. 21: Collision of two solitons with the same parameters 
as in Fig. I2UI but the phase difference of vr between them. 
Solitons bounce at ^ = 180° and 8 = 0°. Finally, the solitons 
merge into a single one, without collapse. 



with a quasi-lD lattice in Ref. Nevertheless, the 

solitons may avoid collapse and separate, if the collision 
time is small due to a large relative velocity, as shown in 
Fig. El. 




FIG. 22: A sequence of snapshots for the collision between 
solitons with individual norms N = 2tt, moving at velocities 
V = ±1 in adjacent troughs (at ro = 5n and ro = 4n) of the 
radial lattice with e = 2. The solitons survive the collision, 
and remain localized, although the collision is inelastic, giving 
rise to radiation loss. 

On the other hand, if the solitons are loosely bound 
(i.e. their norms are small and/or the radial lattice is 
weak), the slow lateral collision gives rise to merger of 
the solitons into one, excessive norm being shed off with 
linear waves, as demonstrated in Fig. 1231 




FIG. 23: The same as in Fig. 1151 but for solitons with smaller 
norms, N = 1.7-k, launched to move at velocities v = ±0.2. 



The existence of such objects was demonstrated by 
means of different variants of the variational approxima- 
tion, and in direct numerical simulations. We have in- 
vestigated the localized states trapped in the central po- 
tential well and in remote circular potential troughs. In 
both cases, a new species, namely, stable radial gap soli- 
tons, have been identified in the model with self-repulsion 
(in addition, the radial gap soliton trapped in a circular 
trough can carry stable pairs of dark solitons on its crest). 
In remote troughs, we have investigated ring-shaped pat- 
terns delocalized in the azimuthal direction, as well as 
strongly localized azimuthal solitons (including moving 
ones). In the attractive model, solutions of both types 
are described by an effective ID equation supplemented 
by the periodic boundary conditions (it is a nonpolyno- 
mial NLS equation with the coordinate running along the 
circumference). In particular, a finite threshold of the 
azimuthal modulational instability (MI) of axisymmetric 
ring-shaped states was predicted, and exact cnoidal-wave 
solutions, generated by the Ml-induced bifurcation were 
found in the framework of the ID equation. Azimuthal 
solitons were also found as solutions to the same equa- 
tion. The existence of stable azimuthally uniform and 
weakly modulated ring-shaped states was corroborated 
by direct simulations. 



Dynamics of completely localized solitons circulating 
in the annular potential troughs was investigated in the 
attractive model by means of direct simulations. In par- 
ticular, the solitons with sufficiently small velocities re- 
main stable indefinitely long, while high velocities give 
rise to leakage of matter into the adjacent (more remote) 
trough under the action of the centrifugal force, which 
eventually destroys the soliton. Collisions between soli- 
tary waves running in the same or adjacent circular chan- 
nels were investigated too. Head-on collisions of the in- 
phase solitons in one trough lead to the collapse; 7r-out of 
phase solitons bounce from each other many times, but 
gradually loose their mutual coherence, and eventually 
merge into a single soliton without collapsing, shedding 
off excess norm. In-phase solitons colliding in adjacent 
channels may also merge into a single soliton. 
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